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Origin of the spectrum of tt electrons that resuhs from coupHng of N Af-\ong acenes via C-C covalent bonding have 
been traced with the use of the Green function formahsm. Exact expressions of acene and graphene Green's functions, 
which are useful for analysis of electronic properties of these macromolecules, are obtained and advanced to the form 
suitable for instructive applications. 
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1 Introduction. Free standing graphene, a monoatomic 
layer of graphite is a new material yji2J that exhibits excep- 
tional electronic quality |[Tl|2][3]|4l|5| and a unique nature of 
charge carriers. On one hand, electrons and holes appear 
as quasi-particles in condensed-matter physics when the 
Schrodinger equation directs our prediction of their behav- 
ior \^J.^dMlE3ilMlMlMM- On the other hand, 
they behave as relativistic particles which are subject to the 
Dirac equation lfT7l[T8l[T9ll20ll2ni22l . Much theoretical ef- 
fort has been focused primarily on the electronic states of 
infinite graphene. In the continuum approximation, these 
states can be mapped into the Hamiltonian of 2+1 dimen- 
sional quantum electrodynamics with Dirac fermions. This 
description is based on the cone-like spectrum near the 
neutrality (zero-energy) points of 2D graphite band struc- 
ture, see II22I and references therein. Less attention has 
been paid to graphene as a molecular species. From this 
point of view, graphene can be considered as a macro- 
molecule obtained by C-C bonding of acene chains as shown 
in Fig. 1. In certain aspects, the macromolecule model pro- 
vides a better approximation for the description of real grapht 
structures |(T6|, particularly, in the nanometer scale, where 
the continuous approximation is of a limited value 1231 . 

As is clear from Fig. 1, Af-\ong linear acenes are ba- 
sic units of NxJV sheet of graphene honeycomb lattice, 
whereas Af-\ong cyclacenes play the same role for zigzag 
{J\f,Q) carbon nanotubes (CNTs). In this sense, the graphene 
and acene spectra are intimately interrelated. The electronic 




Figure 1 A/^-long linear acene as a basic unit of NxJ\f 
sheet of graphene honeycomb lattice (up) and J\f-long cy- 
clacene playing the same role for zigzag (Affi) carbon nan- 
otubes (CNTs). 

structure of acenes, the semi-empirical description of which 
goes back to the Pariser and Parr classic works II24II25L 
is understood fairly well ||26l|27||28l . This study elucidates 
the relationship between the graphene and acene spectra. 
Keeping in mind subsequent applications, our analysis has 
been performed with the use of the Green's function method. 
All derivations are based on the standard operator equation 

G"^ =[EI -H"y\ (1) 
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where is the tight-binding Hamiltonian operator of 
acene (M=A) or graphene (M=G); the electron energy E 
is in unites of \t\, t denotes C-C hopping integral. For the 
Green's function matrix elements, notations {ma\G-^\m' (3) 
^ Gia,m'p and {mna\G°\m'n'P) ee GG_,„,„,^ will 
be used for acene and graphene, respectively; \ma) and 
\nina) denote 2pz orbitals at the ath carbon atom of the 
mth and mrt th hexagon, respectively, as shown in Fig. 1. 
The formal description of acenes and cyclacenes, and gra- 
phene and CNT is very similar. However, in the interests 
of simplicity, only acene-graphene pair will be in focus. 
An extended discussion that also includes cyclacenes and 
zigzag CNTs can be found in Ref. 1301 

2 Acene Green's function. For linear acenes, the ma- 
trix form of solution to Eq. ([T]i is given by ||29ll 



G; 



ma,m'/3 



j;j^^9i,l3^M^jm)sm{^jm'), (2) 



where = 7rj7(7V+l), j = 1, 2, ...,7V, and 



P 



E [i;2-4cos2(^j/2)-l] ,a=^(3^l,r, 



4cos^(Cj/2) 



a = lir),P = ril), 



(3) 

where Vj = [E^ - 4cos2(^j /2)]2 - E^. For the rest of 
matrix elements, we obtain 



{^} E'^-1 



x^g^^^cos2(0/2)sin[0(m-l/2)]sin[e,(m'-l/2)], 

(4) 

( E[E'^-Acos^{^j/2)+ll a = l3^\,p, 



2£;2-4cos2(e,/2), a = A(p),/3 = p(A), 



(5) 



and 



A/" 



7 Xlfa,/3 cos(Cj/2) sin[^j(m-l/2)] sm(ijm') 



(6) 



, _ J £;2-4cos2(e,/2),a = A(p),/3 = /(r), 



E. 



a = A(p),/3 = r(0. 



. (7) 

As shown in |[T6l . defined in Eq. ^ quantities g^^ = 
gl J. and g] ^ = g], ; determine the dispersion relation for an 
iVxA/'sheet of graphene E{Kj, via an equation 



2gi,,cosK, = l-{gjf + igi,.f 



(8) 



where Kj is the second, complementing S^j quantum num- 
ber 2N values of k^^, u = 0,1, N — 1, for each j clas- 
sify the graphene spectrum into 4AfN j, v levels; see be- 
low. 



E 




Figure 2 Band structure of polyacene, E{K=k~^) (red) and 
E{K=k~) (blue), calculated from Eq. ^ for infinite linear 
acene (left panel) and cyclacene (right panel). Continuation 
to imaginary values of wave vectors, fc^ ±7r ± i(5, < 
5 < 21n((l + \/T7)/4) w 0.5, is shown by dashed lines. 
Solid lines of the right panel for k < tt repeats similar 
calculations of Ref. f2T|. Triangles show energy levels for 
Af=l according to Eq. ( fTTT ). 

By performing summation over j, it is possible to ex- 
press the energy dependence of acene Green's function in 
terms of (dimensionless) wave vectors fc+ and fc^, sub- 
jected to equation 

4cos2(fc±/2) = E{E±l). 



(9) 



For example, G^^^^, (^^^j^^^^, = Gfc+ + Gfc-, where 



G 



Gfc 
sin 



fc± 



2{E±1) sin(A/'+ l)fc±' 



(10) 



Other matrix elements have similar expressions. As seen 
from Fig. 2, for the given energy, either both fc+ and k~ 
are real, or one is real, whereas the other one is imagi- 
nary. In the latter case, matrix elements which refer to the 
opposite sites of the chain have two typical terms: One 

is exponential, ~ exp ^/\E\J\f^ , and the second term 
has a singular character. Hence, the probability of electron 
transmission through acenes (~ '-^]^(a/'+i) ED) can occur 
via partly coherent and partly tunneling mechanisms. Tun- 
neling governed by a single exponential factor, as it takes 
place in conjugated oligomers II31II32I . is not possible. This 
is a reflection of a metallic nature of poly acenes t26i271 . 

2.1 Acene electron spectrum. The acene electronic 
structure is determined by equation — = 0; zeros 
of Vj give a four-band spectrum 



E 



c{v) 



.Hi 



Tl±./l + 16cos2(ej-/2) 



(11) 



Copyright line will be provided by the publisher 



pss header will be provided by the publisher 



3 




E 



Figure 3 Green's function matrix element squared, ^j^^, N=2Q, for semiconducting (blue, green) and metallic (red) 
graphene ribbons; the correspondance between colors and Af values is the same, as in Fig. 4. An increase/decrease of N 
results in a shift of the deep and peak in the left and right panels, respectively, towards smaller/larger energies and also, 
in a larger/smaller number of peaks in the mid panel. Differences between the exact Green's function expansion and its 
approximation by the j*th term (mid and right panels), and by a sum of the j*th and (j*+l)th terms of the expansion (left 
panel) are vanishingly small. Dashed blue line shows an approximation by the j * term. Notations = ^g/^ = ^s/^ = 
7r/[2\/3(A/'+l)] are explained in the text. 



where each of two conduction (c) and two valence (v) bands 
has J\f levels. Note that two extra levels = ±1 do not ap- 
pear in Eq. Q as poles of the acene Green's function. The 
states with these energies correspond to zero wave function 
amplitudes at I and r sites and thus, break the acene chain 
into X-p "isolated" pairs. 

For long chains Af>>l and \E\ << 1, spectrum ( fTTT i 
is well approximated hy E ~ ±0^ with q^^— irfi/JV, fi = 
1,2,... . In the limit Af^oo, the quantum number can be 
considered as a continuous variable, k, and (7,, 

q = TT—k acquires the meaning of k separation from the 
zero-energy point at k = tt. The spectrum of polyacene 
(an infinite acene chain) that comes out from the spectra of 
a linear acene and cyclacene, is shown in Fig. 2. 

3 Graphene Green's function. Similarly to Eq. 
matrix elements of graphene Green's function can be rep- 
resented in the form of an expansion 



"Y '^na,n',3 sini^jm) sm{^jm') 



(12) 

where G^,„^„,^ = 0^+ ,^,^ + GL,„'/3- Here, we restrict 
ourselves by representative examples of matrix elements 
refering to the I and r sites, namely. 



sinK± 



sin{K,fN) - gl^sm[K,f{N ~ 1)] 



, (13) 



where labels + and — correspond to, respectively, "plus" 
and "minus" branches of dispersion relation ([8]l. This rela- 
tion can be rewritten as 



E- 



(v) 



+ {-)^Jl + 4cos2(^j-/2) ± 4| cos(^j-/2) cos(kj-/2)|. 

(14) 

Thus, quantum numbers are determined by poles of 
G^^ n'/3 a function of k^, i.e., by the roots of equation 
sin nfN - gj ^ sin Kf{N -1) = 0, or Ull 

sinnfN = =f2 cos(fj/2) siiiK=^(iV + 1/2), (15) 

which together with Eq. (fl4] i gives energies of ANN j, v 
levels. Additionally, the graphene spectrum contains two 
A^-fold degenerate levels with £' = ±1. 

In most cases, it is an interval near the Fermi energy, 
where interesting physics occurs. Whithin this interval and 
for N,J\f >> 1, we have 



(16) 



where = v/37r/[2(7VH-l)], p, « Af, and the upper 
and lower lines refer to an integer and rational value of 
2(A/'+l)/3, respectively. In the case of integer values of 
2(A/'+l)/3=j*, say, fox M=N* , infinitely long armchair 
graphene ribbons (GRs) have a metallic spectrum. In con- 
trast, for M=N*±1, the lowest energy mode is j*=2{M* 
±l)/3, and the spectrum has a band gap {Eg = 2A^/i), 
as defined by the lower line in Eq. ( fTSI l. 

The Green's functions of GRs with the length N and 
N = Af*, Af*+l and Af* — l have a pronouncedly differ- 
ent appearance; see Fig. 3. In what follows, we concen- 
trate on energy intervals below (above) the bottom (top) 
of the second lowest band of metallic GRs, |£'|<Z\(3, and 
whithin the band gap of semiconducting GRs, \E\<Eg/2. 
For these energies, the Green's functions shown in the mid 
and right panel in Fig. 3 are accurately reproduced by a sin- 
gle member of expansion (fT2] |. namely, by the 7*th term. To 
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Figure 4 Solutions to Eq. dTSI ), shown by intersections of 
black, red, and blue curves, corresponding to N =50, 20, 
and 6, respectively, with three horizontal lines 2 cos^j. /2 
for TV = TV* (red), TV* -1 (blue), and A/'* + l (green); TV* =8. 
Green curve, representing A^=3 (<TV/2) does not have in- 
tersections in the region of imaginary k = iS; see text. 



reproduce the energy dependence of the Green's function 
illustrated in the left panel, two terms are needed, partial 
Green's functions j* and 

As immediately follows from Eq. (fT6] l. the E-k rela- 
tion can be satisfied only by real and only by imaginary 
{n=i6) values of the wave vector for metallic and semi- 
conducting GRs, respectively. As a result, an exponential 
factor 

G(i,r,r - exp (^-2N^{Eg/2)^-E^^ , (17) 

appears in the Green's function of semiconducting GRs. 
This factor governs the probability of single-mode elec- 
tron/hole transmission through a potential barrier; see Ref. 
1341 of this issue. For metallic GRs, the number of poles 
depends on TV as, approximately, [^/SN/ (TV+1)+1/2]. For 
semiconducting GRs with TV=TV* — 1, there are no poles at 
all, but in the case of TV=TV* + 1, there is a single pole under 
the condition A^>0.9V3(TV+l)7r— 1/2 or, approximately 
N>Af/2. These peculiarities are illustrated in Fig. 3. 

To conclude this report, the spectra and Green's func- 
tions of graphene and its building blocks, acenes have been 
discussed in parallel that provides a deeper insight into the 
origin and particularities of graphene electronic structure. 
An exact analytical expression of graphene Green's func- 
tion is found in terms of an expansion over partial Green's 
functions. The character of the Green function singulari- 
ties near the Fermi energy has been examined and shown 
to be qualitatively different for metallic and semiconduct- 
ing graphene ribbons. 

This work was partly supported by Visby program of 
the Swedish Institute (SI). 
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